Renormalized 6 model for quantum phase transitions in systems of itinerant fermions 



O 

o 

(N 
5h . 

< 
(N 



I 

c 

o 
o 



> 

oo 
oo 
m 
o 



00 
O 

> 

X 



P. Jakubczylfl 

Max-Planck-Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany and 
Institute for Theoretical Physics, Warsaw University, Hoia 69, 00-681 Warsaw, Poland 

(Dated: April 2, 2009) 

We study the impact of quantum and thermal fluctuations on properties of quantum phase transitions occurring 
in systems of itinerant fermions with main focus on the order of these transitions. Our approach is based on a set 
of flow equations derived within the functional renormalization group framework, in which the order parameter 
is retained as the only degree of freedom, and where the effective potential is parametrized with a <p b form 
allowing for both first and second order scenarios. We find a tendency to turn the first order transitions within 
the bare model into second order transitions upon accounting for the order parameter fluctuations. We compute 
the first and second order phase boundary lines T c as a function of a non-thermal control parameter a 2 in the 
vicinity of a quantum phase transition. We analyze crossovers of the shift exponent if/ governing the shape of the 
T c line when the system is tuned close to a quantum tricritical scenario, where a second order phase transition 
line terminates at a quantum tricritical point. 

PACS numbers: 05.10.Cc, 73.43.Nq, 71.27.+a 



I. INTRODUCTION 

Quantum phase transitions in systems of itinerant fermions 
continue to attract considerable interest * 'i 2 ' 3 ' 4,5 On one hand 
this attention is due to the relevance of quantum critical points 
for understanding collective phenomena like high-r c super- 
conductivity and other non-conventional properties of experi- 
mentally investigated compounds. 6 On the other, it is related 
to peculiarities that so far prohibited the emergence of a fully 
convincing theoretical description of most systems exhibiting 
quantum criticality^ 

The standard approach in describing quantum critical phe- 
nomena in systems of itinerant fermions - the Hertz-Millis 
theory 7 8 - relies on an order parameter field, which is intro- 
duced by a Stratonovich-Hubbard transformation. In a sub- 
sequent step the fermionic degrees of freedom are integrated 
out. The resulting representation of the partition function as a 
path integral over a bosonic field is usually useful only pro- 
vided the action may be expanded in the order parameter, 
and the relevant vertex functions can be evaluated for the val- 
ues of momentum and frequency that correspond to the an- 
ticipated instability. The validity of such an expansion was 
questioned for both zero and finite momentum instabilities in 
the magnetic channel. In case of ferromagnetic transitions it 
was argued that the correct effective action involves additional 
terms 91011 leading to a first order transition, 12 while for the 
case of antiferromagnetic transition in d — 2 integrating out 
gapless fermionic modes can yield singular vertex functions 
all of which were argued to be marginal. 13 Despite these theo- 
retical insufficiencies, Hertz-Millis theory is successful in ex- 
plaining a number of nonconventional properties of many sys- 
tems exhibiting quantum criticality. 6 In recent years this the- 
ory was further extended to account for a number of systems 
not described by the original approach. These include metam- 
agnetic quantum critical points, 14 field-tuned quantum critical 
points where a term describing precession of the order param- 
eter has to be retained in the action, 15 and phase transitions 
induced by a nonequilibrium drive 

An interesting issue concerns the actual order of specific 



quantum phase transitions. Hertz-Millis theory by construc- 
tion allows for continuous transitions only. However, micro- 
scopic models, which usually rely on fermionic degrees of 
freedom and meanfield-like treatment, often predict a first or- 
der scenario in case of both magnetic and charge instabilities 
occurring at Q = 0, and for specific cases corresponding to 
Q + 0^2 Theoretical studies of the interplay between first and 
second order scenarios are further motivated by experimen- 
tal results, where the order of the quantum phase transitions 
in specific compounds is altered by varying external magnetic 
field /l 1819 Upon increasing h, the tricritical point separating 
the first order transition at T < T'" from the second order 
at T > T tn , is shifted towards lower values of T and finally 
vanishes at T — 0. Therefore, the order of the transition at 
small T can be controlled by varying h. Of particular inter- 
est is the case T m = 0, where the transition is second order 
for all T > 0, but scaling behavior of physical quantities at 
low temperatures is in variance with that predicted by Hertz- 
Millis theory ! 20 ' 21 ' 22 An interesting scenario of quantum crit- 
icality occurs also in systems exhibiting a metamagnetic first 
order transition terminating with a critical end point. 14 ^2?. By 
tuning pressure the critical end point is suppressed to T = 
resulting in a distinct type of a quantum critical point. 

In this work we investigate the possibility of altering the or- 
der of quantum phase transitions by order parameter fluctua- 
tions. Subsequently, we study shapes of phase boundaries, fo- 
cusing mainly on the cases where such change occurs. We rely 
on the conventional bosonic Hertz action adapted to phases 
with broken symmetry and retaining also a <p b term in the ef- 
fective potential. We analyze two- and three- dimensional sys- 
tems where discrete symmetry-breaking occurs at Q = (the 
dynamical exponent z = 3). Our results should apply when- 
ever an action of this type can be constructed. Possible ex- 
amples are magnetic transitions with Ising-like symmetry and 
the so-called Pomeranchuk instabilities 24 25 where the discrete 
point-group lattice symmetry of the Fermi surface is broken. 
For these transitions mean-field studies typically predict a first 
order scenario at T = but second order for slightly higher 
temperatures (for T > r'") ' 26 ' 27,28 We argue that the tricriti- 
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cal temperature is lowered by including the order-parameter 
fluctuations. Under appropriate conditions the tricritical point 
may be suppressed to T — 0, in which case a quantum critical 
or tricritical point is realized. 

Our study is based on the one-particle irreducible version of 
the functional renormalization group (RG) i 29 i 30 ' 31,32,33 ' 34 This 
framework was applied extensively in the context of clas- 
sical critical phenomena, where it yields a unified descrip- 
tion of (9(A0-symmetric scalar models including systems in 
d = 2.— It also provides a suitable framework to treat quantum 
criticality. 35,36 We apply the derivative expansion to quadratic 
order, and the effective potential is parametrized in a form 
allowing for the occurrence of both - first and second order 
transition. An advantage of the functional RG is that it allows 
control over effects occurring at different energy scales, and 
treatment of quantities that need not be universal. 

The outline of this work is as follows: In Sec. II we intro- 
duce the bare action to be applied in the RG procedure and 
provide a phase diagram of the unrenormalized theory. In 
Sec. Ill we discuss the applied truncation of the functional 
RG flow equation and derive equations governing the evolu- 
tion of quantities parametrizing the effective potential and the 
inverse propagator upon reducing the cutoff scale. Sec. IV 
is devoted to results obtained for the renormalized phase di- 
agrams in the case T — 0. In Sec. V we consider T > 0. 
We discuss the computed phase diagrams for cases where the 
phase transition is second order for T > and terminates at 
a quantum critical or tricritical point at T — 0. We also ana- 
lyze the case where a tricritical point situated at T > occurs. 
For d = 3 and a choice of system parameters corresponding to 
proximity to the quantum tricritical scenario we compute the 
crossovers of the shift exponent ft, characterizing the shape 
of the critical line, from the Hertz-Millis value ft HM = 3/4 to 
ft*™ = 3/8. This latter value of the shift exponent is specific to 
quantum tricritical behavior in d = z = 3. The obtained value 
of ft 1 " is justified by invoking a scaling argument in Sec. VI, 
where we also classify and discuss scaling regimes emergent 
from our analysis. Finally, Sec. VII contains a summary of 
the work. 



II. BARE ACTION 

In this study we rely on the conventional description of 
quantum critical points in systems of itinerant fermions in 
terms of the bosonic Hertz action* 7 . This is derived by ap- 
plying a Stratonovich-Hubbard transformation to the path- 
integral representation of the partition function of a fermionic 
system. 37 This way the order-parameter field is introduced and 
the original fermionic degrees of freedom are decoupled al- 
lowing one to integrate them out. Subsequently, one expands 
the resulting action in powers of the order-parameter field 
keeping only the most relevant dependence of the propagator 
on frequency and momentum, and neglecting this dependence 
in case of higher-order vertices. In the case of instabilities oc- 
curring at wavevector Q = 0, which is considered here, this 



leads to the following action 

Here <p is the scalar order parameter field and <p p with p = 
(p, oj„) its momentum representation; a>„ = 2nnT with integer 
n denotes the (bosonic) Matsubara frequencies. Momentum 
and energy units are chosen such that the prefactors in front of 
and p 2 are equal to unity. The action is regularized in the 

ultraviolet by restricting momenta to |p| < Ao. The term 
effectively accounts for overdamping of the order-parameter 
fluctuations by fermionic excitations across the Fermi surface. 
The expression Eq. ([TJ is valid for |p|, and ^ sufficiently 
small, which is the limit relevant for the physical situation 
considered here£ 

The potential U[<p] is usually parametrized by a quartic 
form with a single minimum at = or with two minima 
at — +</>o + 0. As we wish to allow for both, first and 
second order transitions, we expand t/[0] to sixth order in 0, 
imposing 

t/[0] = J dT J d d x [a 6 (f> b + a 4 (f) 4 + a 2 (f> 2 ] . (2) 

We consider > to stabilize the system at large \<p\. 

In Eq. (O we assumed the existence of an expansion of 
t/[0], which correctly captures the properties of U[<p] also 
around the non-zero minima. One should be aware, that in 
many physical situations such expansion necessarily has to 
involve also terms of higher order in 0, or simply cannot be 
constructed. The latter, more severe case occurs whenever 
the radius of convergence around = is smaller than the 
distance between = and the minimum at = 0o + 0. 
Our belief is that except for the vicinity of the tricritical point, 
such complication is rather generic. It arises for example in 
case of models describing symmetry-breaking Fermi surface 
deformations in d = 2 (see e.g. Ref. 1281) . While the model 
defined by Eq. © is one of the simplest correctly accounting 
for universal critical and tricritical properties, its applicabil- 
ity in the context of specific microscopic models has to be 
verified for each individual case in the regime of parameters 
where a first-order transition occurs. Also note, that when U is 
expanded around the non-zero minimum 0o, terms involving 
odd powers of field are generated. The coefficients of these 
terms can be expressed as functions of the coupling a 2 , a\, 
«6. There arise just three independent couplings because of 
the assumed inversion symmetry t/[-0] = U[<p]. This model 
does not apply to physical situations where the symmetry of 
U is explicitly broken by including a field coupling linearly to 
0. For example the present work does not apply to the meta- 
magnetic transitions.^?. 

The bare effective potential is now analyzed as a function 
of the parameters a%, ct\, a.(,. It exhibits a single minimum for 
fl2 > 0, a 2 A < 3a2«6, or a 2 > 0, a\ > 0. For a 2 < there are 
two minima located at 



± 0o = + J— -«4 + Ja 2 4 - 3a 2 a 6 . (3) 
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For «4 < 0, a 2 A > 3a2ae, one encounters three minima at <p = 
and = ±0 O . If in addition a 2 A > 4a2«6 holds, U[(po] < U[0] 
and a configuration corresponding to +<f>o is stable, otherwise 
U[</>o] > U[0] and = is the stable state. Along the line 
fl2 = 0, fl4 > a second-order phase transition occurs, while 
a first-order transition line is located along the parabola a 2 A - 
4a2«6 f°i" fl 4 < 0. The tricritical point occurs at — a\ — 0. 
This simple analysis of the bare phase diagram is summarized 
in Fig.l. 
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FIG. 1: (Color online) Schematic phase diagram of the bare if) 6 
model. The line ci2 = 0, a± > marks the locus of second-order 
phase transitions, while the curve a 2 4 = 4a2«6, «4 < corresponds to 
a first order transition. The tricritical point is located at a 2 = a 4 = 0. 
The curves a 2 = and a 2 A = T>aia b for 114 < mark boundaries of the 
region where metastable states are present. 



In course of the subsequent analysis we investigate the sta- 
bility of this phase diagram upon including fluctuations, de- 
pending on the system dimensionality (for d — 2, 3) and tem- 
perature (for T > 0). 



III. FLOW EQUATIONS 

Our analysis is based on approximate flow equations de- 
rived by truncating the exact functional RG evolution equa- 
tion. We use the one-particle irreducible scheme. The flow 
equation22 



d a 1 d A R A 

dA T M = 2V>M + *a 



(4) 



describes the evolution of the effective action T A [0], the gen- 
erating functional for one-particle irreducible vertex functions 
in the presence of an infrared cutoff at scale A, under reduc- 
ing the cutoff scale. In Eq. (O R A denotes the cutoff function 
added to the inverse propagator to cut off modes with mo- 
mentum below the scale A, and F (2) [<p] = 6 2 T A [<p]/6(f) 2 . In 
momentum representation the trace sums over momenta and 
Matsubara frequencies, that is Tr = T 2 M f 4rh ■ ^ n wnat fol- 
lows, we use a regulator term with the optimized Litim cutoff 



tf A (p) = z(A 2 -p 2 )#(A 2 -p 2 ) , 



(5) 



where Z is the wave function renormalization. 

The effective action interpolates smoothly between the bare 
action Eq. (fl} for A = Aq, and the full effective action for 
A — > 0. The quantity T A [if)] can be interpreted as the Gibbs 
free energy of the rapid modes, i.e. modes with momentum 
above the cutoff scale A. By decreasing A fluctuations of 
lower momentum are included and in the limit A — > the 
full Gibbs free energy functional is recovered. The general 
strategy employed in this paper for computing phase diagrams 
amounts to evaluating limA^oF A [0] as function of the (bare) 
parameters describing the system. The equilibrium order pa- 
rameter is extracted by finding the global minimum of the 
resulting free energy. Eq. (HJl describes the flow of an infi- 
nite number of couplings and to make progress we shall im- 
pose a suitable parametrization in terms of a finite number of 
variables. Namely, we apply a truncation in which the effec- 
tive potential U[if>], that is T A [tp] evaluated for momentum- 
independent field, is taken as a sixth-order polynomial as in 
Eq. ©, with flowing coefficients, while the inverse propaga- 
tor is assumed to maintain the form due to Hertz, as in Eq. ([TJ, 
where the coefficients in front of ^ and p 2 may be renormal- 
ized when necessary. Thus, we put 



G- 1 (p,w„)=r (2) [4> = <to] =z w 



Ipl 



•+Zp 2 +2a 2 +fl A (p) ■ (6) 



The quantities Z, Z u depend on the scale A only. Also note, 
that the allowed forms of the effective potential, as illustrated 
in Fig. [T] are typical to <p 6 like truncations. In particular, if no 
truncation of the functional RG flow equation was performed, 
the effective potential should obey bounds on the curvature at 
p = 0, 30 and become convex for A — » 0. 

An analogous truncation (retaining terms up to quartic or- 
der in the effective potential) was applied in Ref. 36, where 
it yields a simple extension of Hertz-Millis theory to phases 
with broken discrete symmetry, capturing the non-Gaussian 
fixed point behavior at finite temperatures. 



A. Effective potential flow 

Upon evaluating Eq. for a momentum-independent field 
<p we obtain an exact equation governing the flow of the effec- 
tive potential U[<p] (see Ref.l30T) 



d,U[p] = 2v d K d l d { W ) 



where p = |0 2 , and 



4W = ^A^r 



d,R A (p) 



Z w ^+Zp 2 +fi A (p) + ZA 2 w 



(7) 



(8) 



Here v d l = 2 d+l 7r d/2 T(d/2), w = w(p) = ^(U'[p] + 
2pU"[p]), and t = log(A/Ao). For our parametrization, the 
partial differential equation Eq. (0 can be projected onto a set 



of three ordinary differential equations. Following Ref. we 
neglect d^Z in d^R A - The neglected terms would yield cor- 
rections to the flow equations that are linear in r\ and there- 
fore irrelevent except close vicintiy of second order transi- 
tion in T > 0. In the latter case these corrections are small 
as compared to other included terms involving 77 due to the 
presence of additional factors involving the interaction cou- 
plings and mass. We use the condition U'[po] = 0, to write 
= f t U'[p ] = d,U'[p ] + U"[p ]d tP Q. Inserting d t U'[p ] as 
obtained by differentiating Eq. (0 with respect to p at p — p , 
one obtains the flow equation for po 



d tPo = 2v d Z- l A d - 2 [3+2 



6a 6 po 



Wlp= 



6fl6P0 + «4 1 

where the threshold functions l d (w) are defined by 

l d (w) = - ±fa) , 



(9) 



(10) 



a 2 
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— 1 <9w 
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From Eq. (0 we conclude that po decreases under the flow. 

This follows from 6fl6Po + a 4 = -^a 2 - 3fl2«6 > 0, and /j(w) > 
0. Analogously, by evaluating the second and third derivatives 
of Eq. © at po and using U"[p] = 48a 6 p + 8a 4 , U'"\p] = 
48<Z6, which follows from the ansatz, Eq. ©, we project out 
the equations governing the flow of a 4 and a§: 



d,a 4 = 12v d A a 



dm = \dv d K d 



4 , (30a6Po + 3a 4 ) 2 _ ,w ^6 
3 2 Z 2 A 4 [ ZA 2 



8^ (30a 6 po + 3a 4 ) 3 
"3 /3 + ^ 



- 6p d,a 6 , 
(12) 

30«6po + 3a 4 



Z 2 A 4 



(13) 

where the threshold functions are evaluated at w correspond- 
ing to po, that is l d n = l d (w)\ p=Po . The flowing coupling a 2 is 
calculated using Eqs. (l9l fT2l[T3l . and Eq. ©, which yields 



a 2 = -4p (a 4 + 3a 6 Po) 



(14) 



The contributions to the flow of the couplings a%, a 4 , a& are 
illustrated in terms of Feynman diagrams in Fig. [2] 

The threshold functions can be split into classical contribu- 
tions, that is contributions from a> = 0, and the quantum parts. 
The Matsubara sums are then evaluated analytically, yielding 
the following expression for the function l d : 



Ad) , Aq) 



i\ = i a ; ci > + 1 



where 



,d(cl) 



2 T 



d(l+w) 2 



(15) 



(16) 



and 



,d(q) 



2T f 1 
f 2 Jo 



.•!l+-^V§-'| . (17, 




FIG. 2: Feynman diagrams representing the contributions to the flow 
equations for the couplings parametrizing the effective potential. 



The 
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Here we substituted y = p 2 /A 2 and f = 
polygamma function \f/\ originated from evaluating Matsub- 
ara sum of the type 



00 

2 ~C) ; 



(18) 



By taking derivatives with respect to w and using i//'„(C) = 
ftn+i (C), we evaluate the higher order threshold functions. For 
any l d one finds a classical contribution proportional to T, and 
therefore vanishing in the zero temperature limit; and a quan- 
tum contribution that survives the limit T — > 0, however, at 
finite T, vanishes quickly as the cutoff scale is reduced. 

One technical remark is in place here. In the derivation 
sketched above it is assumed that U"[po] + 0. This condi- 
tion is not fulfilled whenever the minimum of U[p] at po ± 
disappears in course of the flow, that is the metastable states 
are swept away by fluctuations. Indeed, by inspecting Eq. (0 
we realize that in such case a singularity in <9,po develops for 
some finite t. While solving the flow equations (see Sec. IV, 
V) we will simply terminate the flow whenever such singu- 
larity is encountered and conclude that the system is in the 
symmetric phase with no metastable states present, that is we 
assume that once the metastable states at po are eliminated by 
fluctuations, they would not appear again. 

At this point it is also worthwhile making a comment con- 
cerning some qualitative properties of Eq. (TTZl i. Its right-hand 
side contains terms proportional to powers of (30fl6Po + 3a 4 ) 
and a single negative term involving a^. In the language of 
Feynman diagrams, this term originates from contracting the 
external legs of the 6-point vertex with a single-scale propaga- 
tor (the last diagram contributing to the flow of a 4 in Fig. El). 
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Once we choose the mean field (bare) tricritical point coordi- 
nates po = 0, 04 = as the initial condition of the evolution 
equations Eqs. d9l fT2l[T3] ). only this term survives, and the sys- 
tem is driven into a region of the phase diagram with a\ > 0, 
where a second-order transition occurs for some value of aj_. 
The same happens if one starts with a negative, though suffi- 
ciently close to zero value of the coupling 04. This tendency 
occurs only in the initial stage of the flow. However, it is un- 
likely to be inverted in later stages of the evolution, because 
d t ci4 vanishes quickly as A — » 0, thus, significant renormal- 
izations occur mainly at initial stages of the flow - see Sec. IV 
and V. 



B. Propagator flow 

Within the truncation proposed here, the inverse propagator 
is parametrized with the flowing Z-factors. In what follows, 
we will in addition neglect the renormalization of Z w , As was 
argued in Ref.|3a this is not expected to influence the results 
for the phase diagram. Renormalization of the factor Z will 
be considered only at T > 0, where non-Gaussian critical be- 
havior occurs. The evolution equation for Z is obtained by 
considering the Laplacian of the inverse propagator 



d t Z= —A p [d t G- 1 (v,a) n = 0)] 



p=0 



(19) 



and plugging in the flow equation for G 1 (p): 



d,G-\p) = (12Oa 6 0o + 24a 4 o ) 2 Tr [d t R A (q)G 2 (q) G(q + p)\ . 

(20) 

The latter is obtained by evaluating the second functional 
derivative of Eq. ©. In our calculations we use the order 
parameter anomalous dimension, related to Z by 



1 = - 



c/logZ 
c/log A 



(21) 



By carrying out the Matsubara sums one arrives at the follow- 
ing expression for tj: 



j] = 



2. 2 

- (12Oa 6 0o + 24 «40o) v d TZ- 3 A d - 6 

- 6y 2 (d - l)f fy 2 (*) - 2y 5/2 (2(8 - d) 
(1 + w -y) + 3[(d - 14)(1 + w) + 8y] + 27(1 + w)) 



Tif, 3 (x)-y 3 [2y + W +l] 2 if, 4 (x) 



(22) 



where x — 1 + T~ l (l + w)y 1 ^ 2 . As before, one identifies a clas- 
sical contribution, originating from u — 0, and the remaining 
quantum contribution (involving T ). 



IV. ZERO TEMPERATURE 

In this section we present results for the renormalized effec- 
tive potential at T — 0. As the critical behavior is governed by 
a Gaussian fixed point in this case, 7 we put rj — 0. The classi- 
cal contributions to the threshold functions vanish, while the 
integrals in their quantum parts can be carried out analytically. 
For lf & one obtains: 



l? Q \w) 



ZA J 



1 



1 



nZ^l+w Ud + 3) - 1 



(23) 



as T — > 0. The remaining threshold functions are evaluated 
by taking consecutive derivatives of Eq. ( f23l with respect to 
w. The flow equations are then solved numerically in d = 2, 
and d = 3 for Z = 1, Z a = 1 , and for different choices of ini- 
tial couplings corresponding to a first order transition at mean- 
field level. Namely, we consider the initial value of a& - 1 , 
«4 negative, and 02 ranging between small negative and small 
positive values. For the numerical computations we set the 
upper cutoff An = 1 . The renormalized phase diagrams are 
constructed by inspecting the form of the renormalized effec- 
tive potential (at A — » 0) for the different choices of the ini- 
tial couplings. The locus of phase transition is identified with 
the line in the (02, 04) plane, where the global minimum of U 
becomes nonzero, which happens continuosly or discontinu- 
ously as function of 02 or a\. The metastable states (corre- 
sponding to local minima of U) persist within a region in the 
(«2, fl 4 ) plane in the vicinity of the first order transition line. 
An example solution of the flow equations in case where the 
effective potential is renormalized to a form corresponding to 
a state in the vicinity of a second-order transition, is provided 
in Fig. [3] In Fig.|4]we present a solution for a situation where 



0.07a. 
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+ + + + + ■ 




-0.05 



FIG. 3: (Color online) Solution to the flow equations Eqs. l(9] 1121 
1 1 3I > for d = 2, T = and the initial couplings chosen as follows: 
a 2 (0) = 0.004, a 4 (0) = -0.15, a 6 (0) = 1, and where s = - log A/A . 
The parameter p is renormalized towards lower values, which leads 
to expanding the region of the phase diagram where no symmetry- 
breaking occurs. The coupling 04 increases and reaches a positive 
value (specific to systems exhibiting a second-order transition) in the 
final part of the flow. 
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the absolute minimum of U becomes zero in course of the 
flow, while the metastable states persist for A — > 0. There- 
fore, the average order parameter exhibits a discontinuity as 
function of A. 
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FIG. 4: (Color online) Solution to the flow equations Eqs. ll9l ll2|[T3l l 
for d = 2, T = and the initial couplings chosen as follows: a 2 (0) = 
0.0146, a 4 (0) = -0.25, a 6 (0) = 1, and where s = -logA/A . The 
parameter p corresponding to the global minimum of U exhibits a 
jump as U[po] is shifted above U[p = 0]. 



The results for the renormalized phase diagrams in the vari- 
ables (02, cia) for d — 2 and d — 3 are presented in Figs. [5] |6] 
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FIG. 6: (Color online) Zero temperature phase diagram of the 
renormalized (f> 6 model in d = 3. The renormalized tricritical 
point is shifted towards negative a 4 . The dashed lines and crosses 
mark boundaries of the regions where metastable configurations are 
present in bare and renormalized theory, respectively. 



order of magnitude larger in the case d — 2, where fluctua- 
tions are more pronounced. We also note that the renormal- 
ization of the system parameters occurs at the beginning of 
the flow only. For the present case T = it is caused ex- 
clusively by quantum fluctuations. A consequence of the re- 
sults in Figs. [5] [6] is that there exists a regime of parameters, 
for which the quantum phase transition is first order at mean- 
field level (within the bare theory), and turns continuous when 
fluctuations are accounted for. This change of the quantum 
phase transition's character does not occur in the other direc- 
tion within the framework applied here. 

In the subsequent section, we analyze the system at finite 
temperatures and argue that the tendency towards second or- 
der scenario is even stronger when thermal fluctuations are 
present. Therefore two interesting possibilities are conceiv- 
able - either the transition is second order for T > 0, or it 
becomes of second order for temperatures exceeding a tricrit- 
ical temperature T tn . Of special interest is the case T m = 0, 
where the finite T second order phase transition line termi- 
nates at T — at a quantum tricritical point. 



V. FINITE TEMPERATURES 



FIG. 5: (Color online) Zero temperature phase diagram of the 
renormalized (f> 6 model in d = 2. The renormalized tricritical 
point is shifted towards negative 04. The dashed lines and crosses 
mark boundaries of the regions where metastable configurations are 
present in bare and renormalized theory, respectively. 



As anticipated, the region corresponding to the disordered 
phase in the phase diagram is extended in course of renor- 
malization. One also finds that upon including order param- 
eter fluctuations, the tricritical point is shifted towards nega- 
tive values of a 4 . The magnitude of this shift is nearly one 



In this section we discuss renormalization of the effective 
potential U[<f>] in case of finite temperatures, where both quan- 
tum and thermal fluctuations are present. Conceptually, the 
analysis proceeds along the same lines as for T = 0, however 
the threshold functions now contain both thermal and quan- 
tum parts and the integrals they involve cannot be carried out 
analytically. Moreover, the anomalous dimension 77 has to be 
accounted for, as we encounter non-Gaussian critical behav- 
ior in the vicinity of the transition line. We also note, that 
even slightly off criticality, the presence of anomalous dimen- 
sion influences the results for the phase diagram, in particu- 
lar in d — 2. This happens because rj attains nonzero values 
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at the intermediate stages of the flow, and vanishes only in 
the infrared limit. As a result, the factor Z acquires finite, 
though not negligibly small renormalization. As we checked 
by explicit calculations, a truncation with Z = 1, which breaks 
down only in the immediate vicinity of a second-order transi- 
tion, is not sufficient to determine the order of the transition. 
Namely, at least in d = 2, it gives a very weakly first order 
transition in the regime where the truncation retaining the flow 
of r] predicts a second order scenario. As we also checked, 
the present truncation taking into account the flow of Z repro- 
duces the results of Sec. IV in the limit T — > 0. The relative 
difference in the numerical results for the points exhibited in 
Figs 5,6 using the two different truncations is below 1%. 

Therefore, we analyze the set of flow equations Eqs. (l9l[T2l 
[T3~l l22l with the aim of determining the phase diagrams in the 
variables (02,04) for different temperatures. In a subsequent 
step we reinterpret our results by considering fixed (initial) 
04 < and plotting phase diagrams in the variables (02, T). 
As before, we set the upper cutoff Ao = 1, for the initial con- 
dition we choose = 1, Z = 1, and consider negative initial 
values of 04, and «2 ranging between small negative and small 
positive values. The factor Z w is set equal to unity. An exam- 
ple solution to the flow equations is given in Fig. [7] 
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FIG. 7: (Color online) Solution to the flow equations Eq. {9] [121 
1131 122b for d = 3, T = 0.08 and the initial couplings chosen as 
follows: a 2 (0) = 0.00049, a 4 (0) = -0.05, a 6 (0) = 1, and where 
s = - log A/Ao. This choice of parameters corresponds to the phase 
with broken symmetry, but close to the transition line (within the 
renormalized theory). The coupling a 4 increases and reaches a pos- 
itive value in the final part of the flow. The anomalous dimension 77 
attains nonzero values in the intermediate stages of the flow, although 
the system is slightly off criticality, and vanishes in the infrared limit. 
The system exhibits a second order phase transition upon varying a%, 
although for the same choice of 04(0), (36(0), the transition is of first 
order at T = - see Fig. [6] 



In the solution illustrated in Fig. [7] the choice of 04 cor- 
responds to a situation where, despite the fact that 04 is nega- 
tive, one encounters a second order transition upon varying a^. 
Moreover, for this value of 04, the transition is of first order 
at T = 0, which can be read off from Fig. [6] The initial value 
of fl2 was chosen so that the system is slightly separated from 



the transition, but 77 attains nonzero values at an intermediate 
stage of the flow. 

By performing such an analysis for different values of T, 
«4, and «2, we obtain (for different fixed T), phase diagrams 
analogous to those plotted in Figs|5]|6]for T = 0. Qualitatively 
these phase diagrams are similar, however the position of the 
renormalized tri critical point afp depends on temperature. We 
find that both for d — 2 and d — 3 a'^'(T) is a decreasing 
function, and therefore the tendency towards a second order 
transition is enhanced upon increasing T. In Fig. [8] we plot 
the evaluated a%'(T) for d — 2 and d = 3. The quantity a'£'(T) 
may in general be a complicated function of as well as the 
upper cutoff Ao, and for more sophisticated truncations also 
other parameters characterizing the effective action. There- 
fore, its precise form depends on all these parameters. We be- 
lieve however, that the function would remain monotonously 
decreasing also with better truncations, provided the qualita- 
tive structure of the phase diagram for fixed T is as in Fig. 
1 . This could change for example after including a (p & term 
and allowing for negative «6, which however corresponds to a 
physical situation different to the one studied in this paper. 
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FIG. 8: (Color online) The function a'£'(T) computed from Eqs. l(9] 
El [13] E2l for (a) d = 2 and (b) d = 3. For a 4 > af the system 
exhibits a second order transition upon varying a 2 , while for a 4 < a J' 
the transition is of first order. 



We proceed by evaluating the phase diagrams in the vari- 
ables (0,2, T) for a fixed (initial) coupling a\. For this aim 
one solves the flow equations for different T and a%, inspects 
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the resulting renormalized U and identifies the relevant global 
minimum. We consider two possibilities: a f "(T — 0) < 04 < 
0, in which case the transition is of second order for T > 0; 
and a\ < a'"(T = 0), where the quantum phase transition 
is first order and one encounters a tricritical point at some 
T"' > 0. The results are summarized in Figs [9] \W\ In 
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FIG. 9: (Color online) Transition temperature as function of the con- 
trol parameter a% as computed from Eqs. {9] Q21 Ell I22t in d = 2. 
The ordered phase is situated below the transition line. The plots 
correspond to (a) a 4 = -0.10, and (b) a 4 = -0.21. The transition line 
(a) terminates at a quantum critical point at T = 0, while in case (b) 
the quantum phase transition is of first order, and a tricritical point is 
present at T » 0.006. 



cases where the transition is first order for low temperatures 
we find approximately linear behavior of the phase boundary 
in the limit T — > both in d — 2 and d = 3. However, in 
d = 2 the slope of the transition line is infinite (or very large), 
unlike in d — 3. An important conclusion is that a first order 
transition occurs for T < T' n in the considered situations. We 
do not expect the shapes of the first order transition lines to be 
universal and therefore, in our opinion, it is not quite justified 
to compare the results of the simple (f> 6 model to other calcula- 
tions without first assessing the validity of this approach in the 
considered context. One may note however, that mean field 
calculations of the first order phase boundaries in two dimen- 
sional systems exhibiting discrete symmetry breaking Fermi 
surface deformations yield predictions qualitatively similar to 
our results 28 (very steep, linear behavior of the phase bound- 
ary for T — > 0). Our expectation is that the results of Fig. [9] 
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FIG. 10: (Color online) Transition temperature as function of the 
control parameter ai as computed from Eqs. (l9l ll2l[T3ll22b in d = 3. 
The ordered phase is situated below the transition line. The plots 
correspond to (a) ci4 = -0.03, and (b) 04 = -0.05. The transition line 
(a) terminates at a quantum critical point at T = 0, while in case (b) 
the quantum phase transition is of first order, and a tricritical point is 
present at T ~ 0.05. 



(b) may apply to the system considered in Ref. 28 for nonzero 
uniform fermionic repulsion term u, which is required to as- 
sure that the factor ae in the corresponding Landau expan- 
sion is positive for temperatures reaching down to zero. The 
other requirement is that the tricritical point is located close to 
T — (as in Fig. 4 (b) in the abovementioned work). When 
the results for d = 3 are compared to data from experiments 
on compounds exhibiting ferromagnetic properties (see e.g. 
Refs. I39H40I) . one finds qualitative agreement in that the first 
order transition is located at lower T, reaching down to T — 0, 
while a second order transition occurs at slightly higher T. 
However, another mechanism driving the transition first order 
is present for the case of ferromagnets, and the correspond- 
ing free energy functional involves terms logarithmic in (p. 12 
This fact is not accounted for within the present approach. As 
regards to the situations where a quantum critical point is re- 
alized, we recover approximately linear behavior of the phase 
boundary in d — 2, as dictated by the Hertz-Millis theory. The 
shapes of the phase boundaries in such case for d — 3 are 
discussed in more detail in subsection V B. 
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A. Order parameter exponents 

In this subsection we analyze the behavior of the order 
parameter upon approaching the phase transition line along 
isotherms for which different situations occur. We also com- 
pute the magnitude of the order parameter jump when ap- 
proaching the tri critical point along the coexistence curve and 
extract the corresponding critical exponents. For the purpose 
of this subsection we put a 4 = -0.05 and a 4 = -0.21 for 
d — 3 and d — 2, respectively. This means we focus on the sit- 
uations depicted in Fig.|9](b) and Fig.[l0](b). An investigation 
of the critical exponents characterizing the second-order tran- 
sition in d — 2, 3 in the present context is contained in Ref.l36l 
where the results are also compared to calculations involving 
more sophisticated functional RG truncations 41,42 and exact 
values from the Onsager solution to the Ising model. 

For d = 3 we computed the order parameter, which is given 
by liniA->o0o (<Po being the global minimum of U[<p]) as a 
function of ai, approaching the phase transition line along the 
isotherms T = 0.02, T = 0.04, T = 0.0603 and T = 0.08. 
The results are plotted in Fig.fTTl For T = 0.02 and T = 0.04 
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FIG. 1 1 : (Color online) The order parameter 0o plotted vs a 2 as the 
transition line is approached along the isotherms T = 0.02, T = 
0.04, T = 0.0603 and T = 0.08. For T = 0.02, T = 0.04, where 
the transition is first order, a discontinuity of <f>o occurs at the phase 
boundary. For T = 0.0603, corresponding to the tricritical isotherm, 
the order parameter follows the power law <f>Q ~ (a'" - a^f'", where 
f}, ri as 0.25. Along the isotherm, T = 0.08, one observes a power law 
0o ~ (a" - a 2 y, where [3 = 0.5. Mean-field behavior is observed 
because the plotted points are outside the tiny non-gaussinan regime. 



the critical line. This is not visible in Fig.QT]because the plot- 
ted points are still outside the tiny region where the Ginzburg 
criterion is violated. The crossover to the non-gaussian behav- 
ior as well as sizes of the truly critical region as a function of 
temperature T were studied in detail in Ref . |36| both in d = 2 
and d = 3. 

We proceed by analysing the order parameter jump upon 
varying temperature towards T'" along the coexistence curve. 
The results for d — 2, 3 are plotted in Fig. [12] In d — 3 
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FIG. 12: (Color online) The order parameter jump plotted vs (T' n —T) 
as the tricritical point is approached along the coexistence line for (a) 
d = 3, «4 = -0.05, (b) = 2, «4 = -0.21. In the vicinity of the 
tricritical point the computed order parameter jumps follow a power 
law 0o ~ (T' n - T f. For d = 3 one finds fi » 0.48, being close to the 
mean-field value. For d = 2 /i a 0.079, deviating from mean field 
behavior. 



the order parameter exhibits a discontinuity, as the minimum 
at <p — becomes the global one when the phase boundary is 
reached. The order parameter jump vanishes at the tricritical 
isotherm T as 0.0603, along which the order parameter fol- 
lows the power law <po ~ (a'S ' -aif'", where /?,„■ « 0.25, as an- 
ticipated within mean-field theory. 43 To recover the expected 
logarithmic corrections to the leading power law behavior re- 
quires going very close to the tricritical point, which is not 
achieved here. Along the isotherm, T - 0.08, in the scale ex- 
hibited in Fig.[TT] one observes a power law </>q ~ (a c 2 r - aif, 
where (3 = 0.5 is the mean-field value. In fact, a different be- 
havior, with p 0.31 is expected to occur sufficiently close to 



the jump of the order parameter follows a power law <po ~ 
(T tn - Ty with an exponent fi « 0.48 close to the mean field 
value 0.5. This behavior breaks down at lower temperatures, 
where the system's behavior is no longer expected to be uni- 
versal. For d — 2 a power law with a non mean field exponent 
H as 0.079 is observed. According to our knowledge, the clas- 
sical tricritical behavior in d = 2 was not addressed within the 
functional RG framework. The obtained result compares rela- 
tively well to the accurate value /u ~ 0.094, 43 bearing in mind 
the simplicity of the applied approximation. A better estimate 
of ii would require a more sophisticated truncation of func- 
tional RG. This, and also other aspects of classical tricritical 
behavior in d — 2 remain to be addressed in a separate study. 
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B. Crossover of the shift exponent 

Here we present results regarding the phase boundaries in 
d — 3 in the case, where the transition at T > is sec- 
ond order. We consider two choices of the initial coupling 
a\. a 4 = -0.01, which is separated from a%*(T = 0), and 
a 4 = -0.0302, a value in the proximity of a" A '(J = 0). We per- 
form a careful numerical analysis of the shape of the transition 
lines over a relatively wide range of temperatures and extract 
the effective shift exponent iffff as a function of \ci2 - af^\. 
The quantity ijj e ^ is calculated by fitting a power law to three 
neighboring points of the transition line, which is computed 
numerically. Scaling behavior of the phase boundary shape 
occurs if t]/ e ^ remains constant over a wide range of |«2 -a^l- 
Our findings are summarized in Fig. [13] From Fig.[l3](a) we 
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FIG. 13: (Color online) The effective shift exponent iff" plotted as 
a function of \a 2 - af'\ for (a) a 4 = -0.01, and (b) a 4 = -0.0302. 
Case (b) corresponds to close vicinity of a"'{T = 0). The transition 
line exhibits scaling behavior as dictated by Hertz-Millis theory for 
the case (a). For case (b) a scaling regime with the shift exponent 
t/f ~ 3/8, specific to quantum tricriticality is observed (see the main 
text), and the Hertz-Millis scaling persists only in the narrow vicinity 
of af } , which is not visible in the picture scale. 



read off the standard value of the shift exponent ip = 3/4, as 
dictated by Hertz-Millis theory. However, for case (b) scal- 
ing behavior with ifr a 3/8 is observed. This is related to 
the proximity of the chosen coupling a 4 to the tricritical value 
a%'(T = 0). We also performed analogous computations for 



other values of 04 between -0.01 and -0.03027, yielding the 
observation, that the Hertz-Millis scaling region shrinks as a 4 
approaches a'['(T = 0), and the other scaling regime forms 
sufficiently close to af(T = 0). As a 4 -> a^(T = 0), 
the value of \\i corresponding to the latter scaling platoe ap- 
proaches 3/8. In the subsequent section we provide an ana- 
lytic argument to further clarify the emergent picture and jus- 
tify the obtained value if/' n = 3/8. 



VI. SCALING ANALYSIS 



In this section we show how the obtained values of the shift 
exponents \jj HM and i[>' n can be recovered by invoking a scal- 
ing relation for the free energy in the vicinity of the quantum 
critical point (see e.g. Ref. Qj). The analysis largely follows 
Ref. |3J where if/ HM was derived, and is restricted to d > 2. 
At d = 2 additional logarithmic corrections occur, which are 
not accounted for here. The homogeneity relation for the free 
energy density 

f(a 2 , T, h, a„) = b- d - z f(a 2 b Uv , Tb z , hb y " , a„b M ) (24) 

is written including the magnetic field h and the most rele- 
vant interaction coupling a„. Although the scaling dimension 
[a n ] < 0, it influences the critical behavior because a„ acts as 
a so-called dangerously irrelevant variable^ 4 - Now we differ- 
entiate Eq. ( |24"1 i with respect to h and put h — 0, thus obtaining 
a homogeneity relation for the average order parameter <po at 
h — 0. In the next step we assume, that the dependence of 0o 
on C2 and a n occurs only via the ratio a2/a n . This holds within 
mean field theory, and should not change the value of the shift 
exponent provided this exponent is the same for the transi- 
tion and Ginzburg lines in the vicinity of the quantum critical 
point. This issue was investigated in Ref. 36, where the size of 
the truly critical region was analyzed yielding the conclusion 
that the Ginzburg line provides a very accurate estimate of T c 
in d — 3. Therefore we write 



4> Q (a 2 /a„, T) = b 



-d-z-i 



<>M{a2la n )b llv ~ [a "\ Tb z ) . (25) 



By choosing the scaling factor b such that the first argument 
of 0o on the RHS of Eq. ( 1251 ) becomes a constant, we obtain 



o(a 2 /a n ,T) = (a n /a 2 ) "(f>o(l,T(a„/a 2 ) 



re/(l-v[a„]) 



) , (26) 



with k = v(d + z - - v[a n ]). Now we demand that 

the order parameter </>o vanishes along the transition line 
a 2 (T), that is <po(a 2 (T)/a n ,T) = for any T. This yields 
T{\la 2 {T)) vz ^ l ^ a "^ = const, and therefore the iff exponent 
is identified as 



vz 



1 - v[a„] 



(27) 



Once again we now take advantage of the mean field assump- 
tion identifying the transition and Ginzburg lines and replace 
the exponent v with its mean field value v MFT - 1 /2. Again, 
this step is justified by invoking the results of Ref. 36 and not- 
ing that it suffices to consider a region of the phase diagram off 
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true criticality, in which mean field scaling holds. This yields 



By substituting n — 4 and [04] = 4 — (d + z), we recover the 
formula ft = zl(d + z — 2), derived by Millis^ For d = z — 3 
we obtain ft - 3/4. However, if the most relevant coupling 
corresponds to n — 6, the result is 

if, = z/((2(d + z)-4), (29) 

which for d — z — 3 yields ^ = if/"' 1 = 3/8, as obtained 
numerically in the previous section. 

The corresponding values of the shift exponents can also be 
extracted from Eq. ( 1281 for z = 2, which is not covered by our 
study in Sec. III-V. For this case one obtains ft HM = 2/3 and 
ft ,ri = 1/3. 

Let us also observe, that Eq. (l28l . viewed as a function of 
[a„], generates a sequence of conceivable shift exponents. For 
example for d — z = 3 the sequence is ft n =3/ (2n - 4) for 
n = 4,6,8,.... Monotonous decreasing of this sequence is a 
manifestation of the general fact, that fluctuations tend to sup- 
press the critical temperature. Indeed, as n grows, fluctuation 
effects get less relevant, ft decreases, and the T c line becomes 
steeper. 

Also note, that for d — > 2 + , ft equals 1 for all n and arbitrary 
Z. This follows directly from Eq. (f28). Therefore (up to the 
neglected logarithmic terms), the behavior of the T c line in 
two dimensions is always expected to be linear. 

The analysis performed here and supported by the numeri- 
cal results of Sec. V provides a clear physical picture of the 
system under study. Consider that in addition to the UV cut- 
off Ao there is another scale A'" present. This scale is set by 
(124 - a'")la^. Once A'" > Ao, which always holds for «6 = 0, 
one recovers Hertz-Millis scaling of if/ for |a2 - cif^l ^ Aq. 
If A'" is reduced below Ao, one still observes Hertz-Millis 
behavior as long as the two scales don't become well sepa- 
rated, the only difference being that scaling now occurs for 
|«2 - « 2 0) l Once A'" <k Ao, two scaling regimes are 

present: for \a2 - aL | <k A'" the Hertz-Millis scenario per- 
sists, while for A'" <§: \a% - ai | «C Ao scaling behavior with 
if/'" = 3/8 occurs. If A'" = 0, the quantum critical point is 
replaced with a quantum tricritical point and if/'" determines 
the shape of the whole transition line in the vicinity of T — 
(for| fl2 -4 0) |« A ). 

VII. SUMMARY 

In this work we applied the <p b model to analyze the effect 
of thermal and quantum fluctuations on the phase diagram of a 
system of itinerant fermions exhibiting a quantum phase tran- 
sition at wavevector Q — 0. The analysis is restricted to dis- 
crete symmetry breaking and relies on an effective bosonic 
action analogous to that proposed by Hertz, 7 Renormalization 
of the action is computed from a system of flow equations de- 
rived by truncating the exact functional RG flow equation in 



the one particle irreducible scheme. The flow equations cap- 
ture quantum and thermal fluctuations on equal footing and 
are applicable for T > in d — 2, 3, also in the vicinity of the 
first or second order transition line. 

By analyzing the renormalized phase diagram at T — we 
find that a quantum critical point may be realized even if the 
bare action corresponds to a first order transition. Such sce- 
nario occurs for slightly negative quartic coupling 04. 

We analyzed the dependence of the tricritical quartic cou- 
pling a'" as a function of temperature. The obtained decreas- 
ing form of a'"(T) implies enhancing the tendency towards 
a second order transition upon increasing the temperature. 
Therefore, if the transition is first order at T — 0, it becomes 
second order at higher temperatures, namely for T > T tn . 

We computed phase diagrams in the variables (a2,T), 
where acts as a non-thermal control parameter. In d — 2 
and d — 3 we considered cases where the quantum phase tran- 
sition is first (04 < a'") and second (04 > a'") order. In the 
former situation, the phase boundaries in d — 2, 3 show linear 
behavior in the vicinity of T — 0, where the slope in d — 2 is 
infinite (or very large). Approximately linear behavior of the 
r c -line is also recovered for a\ > a'P in d = 2. 

Subsequently, we analyzed the case 04 > a%' in d = 3. Nu- 
merical computations of the T c line for different choices of 
a\ reveal that the phase boundary follows a power law with 
the shift exponent if/ = 3/4, as predicted by Millis, 8 as long 
as «4 is separated from a'"(T = 0). Upon decreasing 04 to- 
wards a'£'(T = 0), the region in the phase diagram, where this 
scaling is observed shrinks. This is accompanied by forma- 
tion of a different scaling region (at slightly higher T), where 
if/ = if/' n =3/8. When ct\ = a'"(T = 0) the quantum critical 
point is replaced by a quantum tricritical point and the shape 
of the whole transition line in the vicinity of T = is de- 
scribed by if/'". These observations are understood by invok- 
ing scaling analysis relying on mean-field like assumption, by 
which the Ginzburg and transition lines in the vicinity of the 
quantum critical point are identified with each other. By de- 
riving a formula for if/ as function of the dynamical exponent z 
and scaling dimension [a„] of the most relevant coupling, we 
recover iff = 3/4 for quantum criticality in d - z = 3, and 
if/ — 3/8 for quantum tri criticality. 

An interesting avenue for future research would be to re- 
consider these results in the context of specific microscopic 
fermionic models, from which coefficients of the Landau ex- 
pansion can possibly be extracted as functions of physical pa- 
rameters like T and the chemical potential fi for instance. One 
could then check, in what conditions quantum criticality (and 
tricriticality) induced by order parameter fluctuations is con- 
ceivable, i.e. one encounters a%' < a\ < in the vicinity 
of the transition. An interesting candidate for such a study is 
the so called /-model^S. of discrete symmetry-breaking Fermi 
surface deformations, for which z — 3, and which for small 
T and within mean field approximation, exhibits a phase di- 
agram qualitatively similar to the one obtained here in Fig. [9] 
(b)^ 
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